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We introduce a novel kind of spin liquid, the spin 1 chirality liquid, which 
provides a generic paradigm for a disordered spin 1 antiferromagnet in two 
dimensions. It supports spinon and holon excitations, which carry a chirality 
quantum number. These excitations obey half-fermi statistics. The sign of 
the statistical parameter is determined by the chirality. The spin 1 chirality 
liquid is the first example of a two dimensional quantum state which supports 
excitations with fractional statistics but does not violate P or T. 

PACS numbers: 05.30.Pr, 71.10.-w, 75.10.-b, 75.90,+w. 



1. INTRODUCTION 

The key thought of this article is easily expressed. Whenever we have 
a spin liquid (or valence bond solid) state with spin-s on a lattice, we can 
obtain another spin liquid with spin-ms, where m is an integer, by combin- 
ing m copies of such a state and projecting out the spin-ms representation 
contained in 

s s . , , s = ms m (m — l)s 9 ... (1) 



at each site. This is particularly useful if the spin-s state (spontaneously) 
violates one or many discrete symmetries, which we desire not to be violated 
for the spin-ms state we construct. We then only have to combine one copy 
of each of the "degenerate" spin-s states. If there are now m of them, we will 
obtain a state with spin-ms which does not violate any of these symmetries. 

* Dedicated to Peter Wolfle on the occasion of his 60 th birthday. 
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As we shall see below, the projection can be formulated very elegantly using 
Schwinger bosonsB. 

Affleck, Kennedy, Lieb and Tasaki (AKLT)I3 have constructed a series of 
states by availing themselves of this principle for valence bond solids, which 
discretely violate various lattice symmetries depending on lattice type and 
dimension. In this article, I advocate that the principle may have other 
interesting applications, and introduce a particular and generically new ex- 
ample of a two dimensional spin 1 liquid, which I call the spin 1 chirality 
liquid. This liquid supports charge neutral spin-^ spinon (and spinless charge 
1 holon) excitations, which carry a chirality quantum number. This number 
can be x = +1 or — 1- The excitations obey half-fermi statistics, both in 
the sense of Haldane's generalized exclusion principle^ and in the sense of 
the Berry's phases encountered by the state as particles adiabatically encir- 
cle each othera. The chirality quantum number determines the sign of this 
phase. This liquid is the first example of a state with fractional statistics 
excitations in more than one space dimension which does not violate the 
discrete symmetries parity (P) and time reversal (T). 

The spin 1 chirality liquid is, at least at the time of Peter Wolfle's 
60 th birthday, mainly of conceptual importance. As alluded to above, it is 
rather straightforward to formulate it. What is not so easily explained is 
the conceptual importance, the motivation. Why are we interested in spin 
liquids in dimensions larger than one? Why should we believe that they 
generically support spinon excitations, which carry spin ^ and no charge? 
Why do I put forward the hypothesis that spinons in any dimension must 
obey half-fermi statistics? S 

If the reader agrees with my suggested answers to these questions, he 
or she will be quick to realize that spin liquids in two space dimensions 
must either violate P and T (as it has been advocated by R.B. LaughlinQ) 
or support spinon (and holon) excitations, which carry a chirality quantum 
number. (The individual excitations then still violate P and T, but this is 
not any more significant then the P and T violation of right and left moving 
fermions in one dimension: one can always use the complex eigenfunctions of 
any supposed real Hamiltonian H to construct real ones, and hence combine 
them into simultaneous eigenfunctions of H, P, and T.) 

In this article, I present the first spin liquid of this second kind. It is a 
spin liquid with spin s = 1. By now, I have also succeeded in constructing 
a P and T invariant spin liquid for s = s, which I believe to be relevant to 
the problem of CuO-superconductivity£2l. The construction of a generically 
invariant spin liquid is, however, only simple in the s = 1 case discussed here. 
It is one of the more inspired applications of the simple projection method 
outlined above, and at the same time the simplest paradigm for a sophisti- 



S — 1 spin liquids: broken discrete symmetries restored 



cated principle — the principle of fractional quantization and emerging parti- 
cles without P or T violation in two-dimensional antiferromagnets — known 
today. 

In the following three sections, I will illustrate the general projection 
principle with the simplest example, the spin 1 AKLT chain, which is ob- 
tained by combining two dimer or Majumdar-Gosh chainsO. In sections [5] 
and I will introduce the concepts of Jastrow-type wave functions and 
spinon excitations in antiferromagnetically correlated spin liquids in a brief 
review of the Haldane-Shastry modellll, and illustrate why I believe that 
spinons must obey a fractional exclusion principle or fractional statistics, 
regardless of the model specifics includine the dimension. In section [Tj I 
will review the chiral spin liquid stateQ'Ofl, which supports spinon excita- 
tions with half-fermi statistics, but also violates P and T. I will introduce 
the spin 1 chirality liquid, which does not violate any symmetry, in section 
[8] In section |9], I will present a modest amount of numerical work on spin 
1 systems. Finally, I will address the question of spinon confinement in the 



spin 1 chirality liquid in section 10 



2. THE MAJUMDAR-GOSH MODEL 

Majumdar and Gosh0 noticed in 1967 that on a linear spin S = \ chain 
with an even number of sites, the two valence bond solid or dimer states 

l^ G cn ) = n (4 T 4 + u-4i4 + i T )io) = 

(odd) 

v ' i even 

(i odd) 

{I o — o o — o o — o o — o o- ) "even" 
(2) 
| -o o — o o — o o — o o — o ) "odd" 

where the product runs over all even sites i for one state and over all odd sites 
for the other, are exact zero energy ground states of the parent Hamiltonian 

h m G = f Sl S l+1 + ls { S i+2 + I) . (3) 

i 

The proof is exceedingly simple. We rewrite 

HMG = \Y, Hi with ^ = (^ + ^+i + ^+2) 2 -| (4) 

i 

Clearly, any state in which the total spin of three neighboring spins is ^ will 
be annihilated by Hi. (The total spin can only be | or I, as i (g> 4 (8) i = 
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§ © \ ® \-) In the dimer states above, this is always the case as two of 
the three neighboring spins are in a singlet configuration, and (g) \ = \. 
Graphically, we may express this as 



Hi | o— o o ) = Hi | o c^d ) = 0. (5) 

As Hi is positive definite, the two zero energy eigenstates of Hmg are also 
ground states^. 

Is the Majumdar-Gosh or dimer state in the universality class generic to 
one-dimensional spin-i liquids, and hence a useful paradigm to understand, 
say, the nearest-neighbor Heisenberg chain? The answer is clearly no, as the 
dimer states @ violate translational symmetry modulo translations by two 
lattice spacings, while the generic liquid is invariant. A useful paradigm for 
the generic liquid is provided by the Haldane-Shastry model, which we will 
review in section [5] 

Nonetheless, the dimer chain shares some important properties of this 
generic liquid. First, the spinon excitations!^! — here domain walls between 
"even" and "odd" ground states — are free (rather than confined). Second, 
there are (modulo the overall two- fold degeneracy) only M + 1 orbitals avail- 
able for the spinons if 2M spins are condensed into dimers or valence bond 
singlets. This is to say, if there are only a few spinons in a long chain, the 
number of orbitals available to them is roughly half the number of sites. 
This can easily be seen graphically: 



O O O O <j) o — o o — o o — o <p o — o o — 

even odd 



If we start with an even ground state on the left, the spinon to its right must 
occupy an even lattice site and vice versa. The resulting state counting is 
precisely what we will find in the Haldane-Shastry model, where it is directly 
linked to the half-fermi statistics of the spinonsa. 

The dimer chain is further meaningful as a piece of a general paradigm. 
The two degenerate dimer states (]2|) can be combined into an s = 1 chain, 
the AKLT chain, which serves as a generic paradigm for s = 1 chains which 
exhibit the Haldane gapEzUlI. To formulate the AKLT model, we will avail 
ourselves of Schwinger bosons, which we review first. 
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Schwinger bosonaH'Hj constitute a way to formulate spin-s representa- 
tions of an SU(2) algebraEl. The spin operators 

S x + iSy = S+ = a)b 

S x -iSy = S- = tfa (6) 
S z = \{a)a-tfb) 

are given in terms of boson creation and annihilation operators which obey 
the usual commutation relations 

M] = M] = i 

[a, b] = [a,6t] = [at, b] = [at, fet] =0. 

It is readily verified with (0) that 

[S\ S j ] = ie ijk S k where i, j, k = x,y, or z. (8) 

The spin quantum number s is given by half the number of bosons, 

2s = a t a + 6 t 6, (9) 

and the usual spin states (simultaneous eigenstates of S 2 and S z ) are given 
by 

M= {a} p |o). do) 

y [s + my. y (s — m)\ 
In particular, the spin-^ states are given by 

| T ) = a t|0)=c{|0) ||)=6t| ) = c t|o), (11) 

i.e., a) and &t act just like the fermion creation operators cj and cj in this 
case. The difference shows up only when two (or more) creation operators act 
on the same site or orbital. The fermion operators create an antisymmetric 
or singlet configuration (in accordance with the Pauli principle), 

|0,0) = c\c\\0) , (12) 

while the Schwinger bosons create a totally symmetric or triplet (or higher 
spin if we create more than two bosons) configuration, 



|1,1) = ^(at) 2 |0) 




1 1, 0) = aWlO) 


(13) 


1,-1) = ^(6t) 2 |0). 
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Accordingly, may rewrite the Majumdar-Gosh states as 

|«.> = II (44+1 -b\4 +1 ) |0) (14) 



(odd) 
v ' i even 

(i odd) 



(odd) 

This formulation readily suggests a generalization to higher spin. 

4. THE AKLT CHAIN 

Affleck, Kennedy, Lieb and TasakiH noticed that the valance bond solid 
state 

iv> AKLT > = n(«M + i- 6 H + i)io) = 



*-[at,&t] .*2°[at,&t] |0) 
projections onto spin s = 1 



(15) 



is the exact zero energy ground state of the spin 1 extended Heisenberg 
Hamiltonian 



ff AKLT = £ (s^+i + ~ (s>s t+1 ) 2 + fj 



(16) 



with periodic boundary conditions. Each term in the sum (pi]) projects 
onto the subspace in which the total spin of a pair of neighboring sites is 



s = 2. The Hamiltonian (16) thereby lifts all states except (15) to higher 
energies^. The valance bond solid state ( |T5| ) is a generic paradigm as it 
shares all the symmetries, but in particular the Haldane spin gapO'EJ, of the 
spin-1 Heisenberg chain. It even offers a particularly simple understanding 
of this gap, or of the linear confinement force between spinons responsible 
for it, as illustrated by a cartoon: 

— — t — — t — — — 



energy cost cx distance 



We have reviewed the AKLT model here as the simplest and the only 
generally known application of the general projection principle described in 
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the introduction. The Haldane gap illustrates that not only the symmetries, 
but also the physical properties, and in particular the low energy properties, 
may be very different for the combined spin-ms states as for the original 
spin-s state. 



5. THE HALDANE-SHASTRY MODEL 

The Haldane-Shastry modelEBEilil is the most important paradigm for 
a generic spin-^ liquid. Consider a spin-| chain with periodic boundary 
conditions and an even number of sites N on a unit circle embedded in the 
complex plane: 

N sites with spin i on unit circle: 
r]a = e i ^ a with a = 0,1,... ,N-1 




The l/r 2 -Hamiltonian 



where \r/ a — rjpl is the chord distance between the sites a and (3, has the 
exact ground state 

hC>= E € s (^z 2 ,...,z M )sX-...-sf M \u_^_j,) (is) 

{z!,z 2 ,... ,z M } all N spins J. 

where the sum extends over all possible ways to distribute the M = f \ -spin 
coordinates Zi and 

M M 
< s (zi, Z2, . . . , z M ) = (zj - z k f Yl Zj . (19) 

j<k j=l 

This state is real, a spin singlet, and has ground state energy — ^ (iV — J^). 

The proof of solution is rather lengthy 12 H 32 l. We will content ourselves 
here to showing that ( |T9| ) is real, and that it is a spin singlet. As for the 
former, we use 

(zj - z k f = -ZjZ k \zj - z k \ 2 , (20) 
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to write 



M 


M 


M 




1%' ^fcl z j%k 




j<k 


j<k 




M 


M 






\Zj ~ Z k \ 2 J\ G ( Z 3 


) 


j<k 


3=1 






+1 a even 




M 


— 1 a odd. 





C S (^1^2,--- ,Zm) 



... ,., , mi 

j<fc j=l 

where 

= = | i (22) 

The gauge factor G(zj) effects that the Marshall sign criteriaE! is fulfilled. 

We now proof that (19) is a singlet. Since S* ot \tp^ s ) = 0, it suffices to 
show that \tp^ s ) is annihilated by S* t ~ t : 

N 

StoM 3 ) = E 5 « E ^ s (^z 2 ,...z M )st 1 ...st M \{...[) = 

a=l {zi,...z M } 
N 

E J2^o S (Va,Z2,... ,z M )S+...S+ M \i...i) (23) 

{z 2 ,...,z M } a=l 

V * 

=0 

since ip^iVat %2,--- , zm) contains powers 77^, rfc, • • • , and 

N 

^2€ = N5 nfi modN. (24) 

a=l 

The elementary excitations for this model are deconfined and only weak- 
ly interacting spinon excitations, which carry spin-| and no charge. They 
constitute an instance of fractional quantization, similar (both conceptually 
and mathematically) to the fractional quantization of charge in the fractional 
quantum Hall effected. Their fractional quantum number is the spin, which 
takes the value \ in a Hilbert space ( |l8|) made out of spin flips S + , which 
carry spin 1. 

To write down the wave for a j spin spinon localized at site r\ a , consider 
a chain with an odd number of sites N and let M = be the number of 
I or J, spins condensed in the uniform liquid. The spinon is then given by 

M M M 

^al(zi,Z 2 , • • • , Z M ) = JJ (»7 a - Zj) J| (Zj - Z k f J| Zj , (25) 
j=l j<k j=l 
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which we again understand substituted into (|18|). It is easy to verify that 
Stomal = ~ h^af an< ^ * na ^ ^tot^at = wmcn shows that the spinon trans- 
forms as a spinor under rotations. 

The localized spinon ( p5| ) is not an eigenstate of the Haldane-Shastry 
Hamiltonian (|l7|). To obtain exact eigenstates, we construct momentum 
eigenstates according to 

TV 

V>* s j>i,z 2 , ... ,z M )= 2C ^af^l, *2, • • • , *Af). (26) 

where the integer m corresponds to a momentum quantum number. Since 
ip^?(zi,Z2, • • • , zm) contains only powers rj%, rj^, . . . , rj^f and 

N 

2>a< = <5mn m ° d ^ ^ 
a=l 

^ml ( Zl > Z2 ' • • • j^Af) wm vanish unless m = 0,1,... ,M. There are only 
roughly half as many spinon orbitals as sites. Spinons on neighboring sites 
hence cannot be orthogonal. 

To make a correspondence between m and the usual wave number q, we 
translate (|2~6]) counterclockwise by one lattice spacing around the unit circle 



T|C S i > = e i9 |C S i >, (28) 

and find 

M -2m 

q = ttM + 7r . 29 

H 2M + 1 v ' 



The spinon dispersion is given by 



1 



q 2 L v 2 



7T N 2 



mod tt, (30) 



as depicted below. 




allowed momenta 
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The construction ( pqj2q ) can be generalized to many spinon states, which 
weakly attract each other. Then the individual spinon momenta labeled by 
mi (I = 1,2,... , L, where L is the number of spinons) are no longer good 
quantum numbers, as they have a channel to decay into. The construction 
can, however, still be used to count the number of available orbitals for 
spinons, which leads us to the subject of fractional statistics. 



6. FRACTIONAL STATISTICS IN SPIN LIQUIDS 

Consider a Haldane-Shastry chain with an even number of sites N and 
an even number of spinons L. According to (|25|),(p^) and (|27|), the number 



of orbitals available to each spinon is given by M + 1, where M = N ^ L is 
the number of f or j spins in the remaining uniform liquid. (In this rep- 
resentation, the spinon wave functions are symmetric; two or more spinons 
can have the same value for m.) The creation of two spinons reduces the 
number of available orbitals hence by one. They obey half-fermi statistics 
in the sense of Haldane's exclusion principle^. (For fermions, the creation of 
two particles would decrease the number of available orbitals by two, while 
this number would not change for bosons.) 

It is very instructive to use this fractional exclusion principle to count 
the dimension of the Hilbert space spanned by the ground state plus all 
possible spinon configurations! In general, the number of ways to place x 
bosons (as the spinon wave functions are symmetric — the spinons are half- 
fermions formulated in a bosonic representation) into y orbitals is given by 

(s + y-D! (31) 



x\ {y-l)\ 



For a configuration with L spinons, we have x = L and y = 2(M + 1) = 
N — L+2 orbitals, where the factor 2 stems from the internal spin degeneracy 
of each spinon orbital. The total number of states is 



^ L\(N-L + 1)\ ~ [62) 

Li — U 

even 

as it should be for a spin | system with N sites. 

We may interpret the Haldane-Shastry model as a reparameterization 
of a Hilbert space spanned by spin flips (lilf) into a basis which consists of 
the Haldane-Shastry ground state plus all possible spinon excitations. The 
reward for such a reparameterization is that a highly non-trivial Hamil- 
tonian in the original basis may be approximately or exactly diagonal in 
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the new basis. This reparameterization is particularly useful as it usually 
instructs us directly about the quantum numbers of the elementary exci- 
tations. In the fractional quantum Hall effect, we have learned from such 
a reparameterization — or trial wave function — that the quasiparticles carry 
fractional charge. For the spin ^ chain, we can learn that the elementary 
excitations are spinons, which carry spin ^ . 

Of course, it is a condition that the Hilbert space dimensions before 
and after the reparameterization match. It is not very difficult to convince 
oneself that the fractional state counting, and hence the fractional or half 
fermi statistics, is a necessary condition to accomplish this for a spin ^ system 
with spin ^ spinon excitations. This leads me to the following hypothesis: 

Spinons must carry fractional or half-fermion% statistics, regard- 
less of the dimension, in order for the states to count up consis- 
tently. 

In one space dimension, fractional statistics can only be defined through 
a fractional exclusion principle, which does not yield any symmetry viola- 
tions. In two space dimensions, however, we can define statistics alterna- 
tively through the Berry's phase acquired by the state as we adiabatically 
exchange particles by moving them counterclockwise around each otheru. 
This phase is given by e tS , where 9 is the statistical parameter. 



m - IV) 

For 9 = ir we have fermions, for 9 = bosons, and for both 9 = ^ and 
9 = — ^ half- fermions. The choice of sign for 9 is physically meaningful, as 
the allowed values for the relative angular momentum I depend on it: 

Q 

I = — Yin where n integer. (33) 

Clearly either choice of sign for half- fermions violates parity (P) and time 
reversal (T). 

If we now assume that both definitions of statistics match for spinons 
in two dimensional spin liquids, as they do for the quasiparticles in the 
fractionally quantized Hall effect, the hypothesis of half fermi statistics leaves 
us with only two choices. The first is that the spin liquid ground state 
violates P and T spontaneously and thereby fixes the sign for the statistical 
parameter for the spinons. This possibility has been advocated by R.B. 
Laughlin. The chiral spin liquidloi, which we will review in the following 
section, provides us with a paradigm for this situation. 
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The second possibility, which I would like to advocate, is that the 
spinons carry a chirality quantum number, which determines the sign of 
the statistical parameter 9. Then there is no need for the ground state to 
violate any symmetries. Spinons of different chiralities map into each other 
under P and T. The situation is somewhat analogous to fermions in one di- 
mension, where we have right and left moving excitations, which are P and 
T conjugates of each other. I will construct the simplest paradigm for such 
a liquid, the spin 1 chirality liquid, in section |8] 



7. THE CHIRAL SPIN LIQUID 

The chiral spin liquidi'B&i may be viewed as a brute-force generalization 
of the Haldane-Shastry wave function to two space dimensionsEJ. Consider 
a periodic one-dimensional lattice on the real axis of a complex plane, with 
lattice points at integer values: 



N-l N(=0) 



Filled and empty circles represent even and odd integers, with gauge fac- 
tors G(z) = —1 and G(z) = +1, respectively. The Haldane-Shastry wave 
functions d2lh then becomes 



M M 

ipf a (zx,. . . , z M ) = JJ G(zj) Yl sm(^( Zj - z k ) 



j<k 



(34) 



where we took advantage of the fact that now the coordinates Zj are real. 

The chiral spin liquid is obtained by extending the lattice from a circle to 
a cylinder, or from a segment of the real axis to a strip in the two-dimensional 
complex plane: 



o • o • o • o 



ooooooooo 



o o o o o o o 
o • o • o • o 



L»(=0) 
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where G(z) = (— l)( a! + 1 )(#+ 1 ) for lattice site z = x + iy. The wave func- 
tion for a chiral spin liquid on this cylinder is given by (|34|) multiplied by 
an exponential factor, which effects that the density of spin flips is \ for 
—Ly/2 < y < Ly/2, i.e., up to the boundaries of the liquid: 

M M , \2M 

^(zi,... ,z m ) = UG(z j ) l[ s J^( Zj - Zk )\ JJe-^l 2 . (35) 

j=l j<k V x J j=l 

The chiral spin liquid takes a more familiar form if we consider open 
boundary conditions. The wave function for a circular droplet of fluid is 
simply 

M M M 

^(z 1 ,...,z M ) = l[G(z j ) l[( Zj -z k f [Je-I^'l 2 . (36) 
j=l j<k j=l 

Note that the exponential in ( [35]) or (|36| ) corresponds to a (fictitious) mag- 
netic field of strength 2-7r/plaquet. Apart from the gauge factors G(zj), ( |36"| ) 
is formally equivalent to the wave function for a fractionally quantized Hall 
liquid for bosons at filling fraction v = |, which corresponds to half a particle 
(or spin flip) per plaquet. t/^ sl is a spin singlet. 

The chiral spin liquid is in general complex as the z'b are complex. 
There are two "degenerate" states, ■0+ SL an d its complex conjugate , i/'° SL , 
which are P and T conjugates of each other. (I have written "degenerate" 
in quotation marks as no-one has ever been able to construct a local parent 
Hamiltonian for this liquid.) Since the chiral spin liquid violates the discrete 
symmetries P and T, it is not in the same universality class as the generic 
disordered spin \ antiferromagnet in two space dimensions, which I presume 
to be stabilized for some lattice type or upon doping. 

It is rather easy to formulate spinon excitations for the chiral spin liquid. 
In analogy to both the quasiholes in fractionally quantized Hall liquids and 
the spinons ( [25] ) for the Haldane-Shastry model, we write the wave function 
for a spinon localized at r] 

M M M M 

^ l (z 1 ,...,z M ) = l[(Va-z j ) [[ G(^) J\{z3-z k f l[e-^\ 

3=1 j=l j<k j=l 

(37) 

The spinons obey half-fermi statistics, both in the sense of Haldane's exclu- 
sion principle as well as in the sense of the Berry's phases encountered by 
the state as we exchange particles by winding them counterclockwise around 
each other. We can evaluate this phase using the adiabatic transport argu- 



ment of Arovas, Schrieffer, and WilczektJ. We obtain tt/2 for (37) and —ir/2 
for its complex conjugate. 
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The chiral spin liquid is a wonderful paradigm, with only one imperfec- 
tion: the violation of P and T. There is no indication that these symmetries 
are spontaneously broken in any two-dimensional spin system. 

Like in the case of the Majumdar-Gosh state discussed in section |2], the 
chiral spin liquid can be used to construct a symmetry invariant spin liquid 
with s = 1, the spin 1 chirality liquid. For this construction, it is propitious 
to formulate the chiral spin liquid in a basis which constitutes of electron 
creation operators rather than spin flips. To begin with, we use M = N/2 
to rewrite (^) as 

M , M \ 

r + SL (z 1 ,...,z M ) = H(G(z J )e-^ 2 J] {Zj-z k )\. (38) 

3=1 ^ k=l ' 

(Mi) 

Let wi with 1 = 1,2,... ,N — M be the lattice sites not occupied by the Zj's. 
Then the octopus theoremil 

N N 

G{ Zj ) e-f W 11 ( Zj ~Va) = Yl e + ^ M2 ■ const., (39) 

a=l a=l 



where the first product runs over all lattice sites except Zj, implies 

M / M \ M M 

*?-(*,. ..,«)= n n -4^r n n 



3=1 \ 1=1 ■> 1 / 3=1 1=1 

(40) 



Let 



S[z,w] = (0 \ c Zli . . . c ZMl c Wli . . . c WN _ Ml \ 11 .. . ... | ), (41) 

all A 1 " spins J, 

be the sign associated with ordering the z's and io's according to their lattice 
positions. Then we can use 

M M N—M N-M 

- n n ^ ~ w n ^ - ^ ■ 

7<fc 7=1 Z=l Z<m 

(42) 

M N—M 

■ Y[e-^ 2 Y[ e -fl^l 2 = <S[2,w] -const. 
i=i /=i 

to rewrite (|^) as 

V>+ SL (zi,... ,zm) = S[z,w] <f>{zi,... ,zm) <f>(wi,... ,wn-m) (43) 
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where 



M M 

Mtv i _ .12 



<j>{z x ,... ,z M ) = Jlte-Zk) l[e-^ (44) 

j<k j=l 

is simply the wave function for a filled Landau level in a (fictitious) magnetic 
field with flux ^^/plaquet. If we rewrite (36) in terms of ( |44| ) and compare 
it to (|43l), we obtain the lattice particle-hole symmetry 



M 

\[G(zj) <f>(zi,... ,z M ) = S[z,w] (f>{wx, . . . ,w N - M ) ■ const., (45) 
i=i 

which holds for any M. The chiral spin liquid ground state, where M = N/2, 
is according to ( [42} ) simply given by 

\l/l+ Sh ) = ^2 <f>(zi,... ,ZM)<f>(wi,... ,Wm) ■ 

{zi,... ,z M ;W!,... ,w M 

} (46) 

' C ziV C zm1 C w 1 1--- C w m i\°) ' 

where the sum extends over all possible ways to distribute the coordinates 
Zj and wi on mutually distinct lattice sites. It is often convenient to write 
as 



|^ SL > = PcwK> (47) 
where the Gutzwiller projector 

N 

P G w = J] C 1 - c h c n 4i c n) ( 48 ) 

i=l 

eliminates doubly occupied sites and (V'sd) is the Slater determinant wave 
function for the lowest Landau level filled once with M = f t -spin and once 
with M j-spin electrons. 

The spin singlet property of the chiral spin liquid mentioned above is 
now evident. Clearly the Slater determinant I^sd) is a singlet. Since the 
Gutzwiller projector commutes with the spin operator on each site, 

[P GW , Si] = 0, (49) 

i^csl^ mus t kg a s i n gl e t as well. 

The Gutzwiller form (|47]) also allows a rather elegant formulation of 
the spinon excitations. For example, a state with L j-spin spinons at sites 
rji, . . . ,rj L is given by 

\^+ml-VLi) = Fewest ' • • c vl1 \^sj^ L ) > ( 50 ) 
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where N + L = 2M must be an even integer. The equivalence of ( |50| ) to ( |37|) 
is readily seen with (^). This form nicely illustrates the (fractional) spin 
2 1 of the spinon. The electron annihilation operators create inhomogeneities 
in spin and charge before projection. The projector enforces one particle 
per site and hence restores the homogeneity in the charge distribution, but 
commutes with spin. We are left with a neutral object of spin 5 . 

Note that the spinon coordinates rj in either ([37]) or ( |50D do not have 
to coincide with lattice points. This allows us to define a winding path and 
thus a statistical parameter via a Berry's phase encountered upon adiabatic 
interchange. 

Holon excitations, which carry a positive unit charge and no spin, are 
constructed from spinon excitations by annihilation of an electron at the 
spinon site, which must now coincide with a lattice point. For example, a 
chiral spin liquid with L holons at sites rji, . . . ,r]L is given by 

|V , +^io...r ?L o) = c f?iJ. • ' ' C VlI Pgw Crjit • ' ' C f]h\ |V'sd +L ) > (51) 

where N + L = 2M must again be an even integer. 

The chiral spin liquid, when formulated as a Gutzwiller projection of 
filled Landau levels, is readily generalized to different lattice types. 

8. THE SPIN 1 CHIRALITY LIQUID 

The spin 1 chirality liquid is obtained by combining t/^ sl and its com- 
plex conjugate ip^sL v j a a p ro j ec tion to spin 1 on each site, which can be 
accomplished elegantly via Schwinger bosons. Let 

|V"+ SL ) = ^2 0(^1, ■■■ ,z M )(l)(w 1 ,... ,w M ) ■ 

{21,... ,z M ;wi,... ,w M } 

■ ai 1 ...ai M bl 1 ...bl M \0)= (52) 

= ^ SL [at,&t] |0>. 

be a chiral spin liquid ( |46|) written in terms of Schwinger bosons. The spin 1 
chirality liquid is then given by 

|V S1CL ) = V c + SL [a\tf] ^ SL [at,&t] |o>. (53) 

This is the central proposal of this article. 

The spinon or holon excitations for the spin 1 chirality liquid are just 
the spinon or holon excitations of the individual chiral spin liquids. As the 
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spinons for the chiral spin liquid are massive, we expect an energy gap for 
the spin 1 liquid as well. The spinons in the latter carry a chirality quantum 
number x-> which is +1 if they are constructed as excitations of ■0+ SL and —1 
if they are constructed as excitations of ■0° SL . 

Let us first consider the case where we only have spinons of chirality 
+. They span a Hilbert space of dimension 2 N , just as in the case of the 
chiral spin liquid or Haldane-Shastry chain. The spinons obey a fractional 
exclusion principle according to half fermi statistics. The adiabatic trans- 
port argument of Arovas, Schrieffer, and Wilczekil remains unaffected by 
the projection as well. We can evaluate the phase acquired by \E'° SL [cr, 6'] 
and hence |^ S1CL ) as we exchange two + chirality spinons by winding them 
counterclockwise around each other, and obtain tt/2. If we had only spinons 
of chirality — , we would obtain —tt/2. In this sense, spinons of either chiral- 
ity obey half-fermi statistics, with statistical parameter 

= x\ (54) 

Non-trivialities arise, however, as we consider states with spinons of 
both chiralities. If spinons of different chiralities were independent of each 
other, the total Hilbert space dimension of the ground state plus all possible 
spinon excitations would be 2 N -2 N = A N . The Hilbert space dimension 
for a spin 1 system with N sites, however, is only 3 N . The spinon basis is 
consequently overcomplete, and spinons of different chiralities cannot always 
be orthogonal. They have a significant overlap when they are spatially close 
to each other. 

What does this mean? One possibility is that a parameterization of the 
spin 1 Hilbert space in terms of a ground state plus spinons excitations is 
not sensible. This would be the case if there is no energy scale at which 
spinon degrees of freedom offer an adequate description of the system. I do 
not believe in this possibility but cannot rule it out at this stage. 

The other possibility is that there is a way to transform the 4^ non- 
orthogonal orbitals into 3^ orthogonal orbitals, which then provide the de- 
sired Hilbert space basis. Note that a similar situation had occurred in the 
Haldane-Shastry model, when we first constructed spinon excitations in real 
space and ended up with a vastly overcomplete basis. The problem was re- 
solved by a transformation to momentum space. Of course, the possibility 
of such a transformation confronts us with a variety of unresolved questions. 
The most pressing one, in my opinion, concerns the statistics. Since we do 
not expect the transformation to dissolve the concept of spinons, it will be 
linear and hence preserve the fractional exclusion principle. It is not clear, 
however, what happens to chirality quantum numbers, winding phases, and 
relative angular momenta. 



M. Greiter 



Fortunately, the issue of overcompleteness only matters when we have 
spinons of different chiralities nearby. For the low energy dynamics of a spin 
1 chirality liquid another question is more immediate: are the spinons free, 
as they are in spin i chains or in chiral spin liquids, or are they confined, as 
they are in spin 1 chains? Even though I will not be able to give a definite 
answer, I will present a modest amount of evidence for a conjecture in section 
|10.| First, however, I will interpret a few numerical studies on the the spin 
1 chirality liquid ground state in the following section. 



9. FINITE SIZE STUDIES 

In this section, we report on a few numerical studies of the spin 1 chi- 
rality liquid (S1CL) for square, triangular and kagome lattices, as well as 
the chiral spin liquids (CSLs) used to construct them. The calculations are 
performed on finite size systems with periodic boundary conditions (PBC). 
To construct the spin liquid states for this geometry, we have to general- 
ize the wave function (E3|) for a filled Landau level to PBCs, which is in 
essence done by replacing the factors (zj — z^) by odd Jacobi theta functions 

(Zj - z k ,r)B 

The most severe limitation for exact diagonalization studies is the size 
of the Hilbert space. For a spin 1 system, the dimension of the subspace 
with S* ot = is given by 



JV! 



n\n\ (N - 2n)\ 

n=0 v ; 



(55) 



which amounts to 5,196,627 for N=16. This limits the system size to an order 
of 16 sites, which is far too small to address most questions of interest. For 
example, it is well established that all the spin 1 Heisenberg antiferromagnets 
we study possess long range order. This feature, however, is hardly manifest 
in the systems we consider here. The comparisons of the spin liquid trial 
wave functions with exact ground states of Heisenberg models reported in 
Table |l| merely indicate to which extent the energetically relevant nearest- 
neighbor correlations of the models are captured by the liquids. If these 
correlations are sufficiently close, we may conjecture that the spin liquid is 
stabilized if we disorder the system by doping. 

Let us now take a look at Table |]. First of all, the results confirm 
that the P and T violating CSLs are inadequate for all three lattice types. 
The energies always differ by an amount of order 10 % or larger, which is 
enormous in light of the fact that the energy of the classical ground state is 
only about 30 % higher than the energy of the quantum model. Let us now 
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Lattice type 


square 


t t*i o n m l qv 
Llldll^UlclI 








J' = 


,V = \J 






Number of sites N 


16 


1 f\ 


1 9 


Principal displacement L\ 


(4,0) 




f4 0) 






(0,4) 


(2 2\/~3) 


f2 2\/3l 


LL bound, phases <pi,<f2 


7T 


TT 


2-7T 47T 

3 ' 3 


7T 37T 
4 ' 4 




rp 

-^cxact 


-11.2285 


-8.4579 


-8 5555 


4448 


S — 2 


-E'cSL 


-8.2572 


-7.6435 


-7 8251 


-4 9756 

'-t.iJ 1 'J\J 




^^jCSL | ^exact ^ 


.6727 


.8626 


.0000 


.0000 




-^cxact 


-38.3047 


-25.0858 


-25.3699 


-17.6208 






-30.7987 


-24.4782 


-24.1504 


-15.7025 




AE/ -E exact 


19.6 % 


2.42 % 


4.80 % 


10.89 % 


S= 1 


^^SICL |^cxact ^ 


.6019 


.9484 


.7948 


.5527 




-^SlFl/2 


-38.1286 


-24.8957 


-25.0184 


-15.9249 




AE/ -Bexact 


0.46 % 


0.76% 


1.39% 


9.62 % 




^SlFl/2 |^cxact ^ 


.9927 


.9784 


.8962 


.02954 



Table 1. Energy expectation values and overlaps for the chiral spin liquid 
(CSL), the spin 1 chirality liquid (S1CL), and a spin 1 liquid constructed by 
combining two identical copies of the exact spin ^ ground states according to 
the Schwinger boson projection principle discussed in this article (S1F1/2), 
for finite lattices of various types with PBC. The trial states are compared 
to exact ground states of nearest-neighbor Heisenberg models with J = 
1 for all three lattice types, plus a model augmented by a next-nearest- 
neighbor coupling J' on the square lattice. The displacement vectors L\, 
1/2 span the principal region for the PBCs in units of lattice constants. 
The boundary phases ip\ arid if2 are a property of the Landau level wave 
functions adapted to PBCsc3 and projected onto the lattice, which are used 
to construct the chiral spin liquids: (p% and if2 are the phases acquired by 
( Pl 1 .l 2 (zi, ■ ■ ■ ,zm) as one of the coordinates (which must coincide with a 
lattice point) is translated by L\ or L 2 , respectively, plus the phase acquired 
by a unit charge under this translation through the coupling to the vector 
potential generating the magnetic fleld.EO These phases have been chosen 
to minimize the energies of the corresponding spin liquids, which slightly 
depend on it. 
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turn to the S1CL for the square lattice. For a Heisenberg antiferromagnet 
with only a nearest-neighbor coupling J, the S1CL is is way off as well, as 
shown in the first column of Table |]. If we use, however, the S1CL as a trial 
wave function for a Heisenberg model with a next-nearest-neighbor coupling 
J' = J/2, we find an energy difference of 2.42 % and an overlap of 95 %, which 
is reasonable. There is a very realistic chance that the S1CL is realized in 
this antiferromagnet if it is disordered by additional frustration or disorder. 
It may further be a candidate for a disordered state on the triangular lattice 
but appears inadequate for the kagome lattice. 

In the last three lines on the bottom of Table |], I have evaluated energies 
and overlaps of yet another trial wave function for s = 1 antiferromagnets, 
which is simply obtained by combining two identical copies of the exact 
ground states for the corresponding s = \ Heisenberg model 

|V cxact ) = * 1/a [a!,tf] |0) (56) 

by Schwinger boson projection: 

|V slF1/2 > = ^> 1/2 [a\b ] ] m 1/2 [a\tf] |0>. (57) 

I call this trial state "spin 1 from 1/2" (S1F1/2). Comparisons of this state 
with the exact ground states of the spin 1 antiferromagnets reveal informa- 
tion about how different the correlations in the spin 1 and spin ^ models are. 
For the nearest-neighbor Heisenberg model on the square lattice, the energy 
of these states differs only by 0.46%, with an overlap of 99%. I attribute 
this to the fact that both states share the same type of long range order. For 
the triangular lattice, this overlap is already reduced to 90 %. I interpret this 
as an indication that the order is much more pronounced in the triangular 
spin 1 system. For the kagome lattice, the overlap is virtually zero, which 
confirms that spin | and spin 1 models possess very different correlations. 
The spin ^ system is an intrinsically disordered spin liquidc2, while the spin 
1 model on the kagome lattice possesses order. Most likely, there is a disor- 
dered state for the spin 1 kagome system which can be stabilized by doping, 
but it is apparently not obtained by Schwinger boson projection of two iden- 
tical spin ^ liquids. I conjecture that nature avails herself of the possibility 
that the individual spin ^ states used in this construction may violate a 
discrete symmetry, which is not P or T, but symmetry under rotations by 
7r/3 modulo 27r/3. This is, however, merely a speculation. 

10. SPINON CONFINEMENT 

In section W\, we illustrated with a little cartoon how a linear confine- 
ment force between spinons arises for the AKLT chain, even though spinons 
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spinon -.552 -.381 -.478 




Fig. 1. Spin-spin correlations in two spin ^ Haldane-Shastry chains with 
15 sites and one spinon each (a) and a spin 1 Gutzwiller chain (b) con- 
structed via Schwinger boson projection from two Haldane-Shastry chains 
with spinons at different sites. The thickness of the lines reflects the mag- 
nitudes of the antiferromagnetic correlations {SiSi+i) on nearest neighbor 
sites. 

are free for the individual Majumdar-Gosh chains. This example might sug- 
gest that if spinons are confined, this property can always be understood 
through such a simple picture. Let me first demonstrate that this is not the 
case, and then turn to the subtle question of whether spinons are confined 
in the spin 1 chirality liquid. 

Consider a spin 1 Gutzwiller chain (S1GW), that is, a spin 1 chain 
obtained by Schwinger boson projection of two identical Haldane-Shastry 
chains: 

|V S1GW ) = *f[at,&t] *f [at,&t] |0). (58) 

As the individual Haldane-Shastry chains are completely disordered, there 
is no indication of confining forces on the level of the cartoon mentioned 
above. Nonetheless, the S1GW chain does not support free spinons. There 
are several ways to see this. The first is that the S1GW chain is in the same 
universality class as theAKLT or spin 1 Heisenberg chain, which is known to 
possess an energy gapc^l. Since the spinons for the Haldane-Shastry are (in 
contrast to the Majumdar-Gosh chain or the chiral spin liquid) massless, the 
energy gap must stem from spinon-spinon interactions, or confining forces. 
This is in accordance with Haldane's analysis^ of the one-dimensional 
nonlinear a model, which shows that a topological term and hence decon- 
fined spinons exist only for half-integer spin chains. 

The second way to see that spinons in the S1GW chain are confined is 
by direct computation of the confining potential, as shown in Fig. HI. The 
numerically accessible system size is again to small for a full analysis, as we 
must consider periodic boundary conditions^ and are hence limited to a few 
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Fig. 2. Spin-spin correlations in a chiral spin liquid (a) and a spin 1 chirality 
liquid (b) constructed with it for a 16 site square lattice with periodic bound- 
ary conditions in the presence of two spinons at the positions indicated. 



lattice spacings when we pull two spinons apart. Nonetheless, Fig. p] clearly 
shows a dimerization of the spin correlations between the spinons. The 
dimerized chain segment compares unfavorably in energy to the uniformly 
correlated chain far away. I should mention here that when I say spinons 
are confined, I do not necessarily mean that there is a string force between 
them, aa it is the case for spinons in AKLT chains or two-leg Heisenberg 
ladders^. I merely mean that it costs less energy to create another spinon 
pair than to completely separate two spinons. In the S1GW chain, the 
confining force approaches a constant for large spinon separations, as a single 
localized spinon in a Haldane-Shastry chain dimerizes the spin correlations 
over the entire chain. The strength of the dimerization and hence the force 
between widely separated spinons, however, decreases with the system size, 
and vanishes in the thermodynamic limit E3 

Let us now turn to the question of whether there is spinon confinement 
in the spin 1 chirality liquid. First of all, from an analysis of the nonlinear 
a model in two dimensions, there is no reason to expect deconfined spinons 
for either half-integer or integer spins.il The fact that the chiral spin liquid 
supports deconfined spinons may be interpreted as yet another indication 
that it does not describe the universality class of the generic spin \ liquid in 
two dimensions, which is believed to be of central importance to the problem 
of CuO superconductivity. 

Evidence for confining forces in the spin 1 chirality liquid is presented 
in Fig. ^. In a CSL with two spinons (a), the energetically relevant nearest- 
neighbor spin correlations are slightly enhanced near the spinons. This en- 
hancement is due to the fact that the sites around the spinons effectively 
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have fewer neighbors, as there is no need to correlate the links to the spinon 
sites. Combining this + chirality CSL with a — chirality CSL ground state 
(which has spatially uniform correlations) yields a S1CL with two + chiral- 
ity spinons (b). The spin correlations for this state are slightly depressed in 
the immediate vicinity of the spinons, in sharp contrast to what one would 
expect from the correlations in the CSLs. I interpret this depression as 
an indication of confining forces between the spinons. I should emphasize, 
however, that the system size is too small to allow for a conclusion. 

If one believes that there is a confining force between spinons in the 
S1CL, one is led to wonder where it is coming form. In the spin 1 chain, 
we have been able to trace it to a dimerization in the spin correlations 
in the chain segment between the spinons. The spin correlations provided 
information about both the existence and the origin of the force. In the 
S1CL, however, Fig. || only suggests that a force exists, but does not reveal 
its origin. It seems to come from no-where. The origin of the force, I believe, 
only reveals itself if one considers chirality correlations, that is, expectation 
values of the chirality operatorQ 

X = Si ■ (Sj x S k ), (59) 

where i,j, and k are three lattice sites on the same plaquet. This operator 
has been used as an order parameter for the CSL, as it effectively provides 
a local measure of the P and T violation in a spin liquid. The expectation 
value is trivially zero in the S1CL ground state. As the individual spinons 
in the S1CL carry chirality and hence violate P and T, we expect a cloud of 
chirality, that is, non-zero expectation values of ( |59[ ) on triangles, around an 
isolated spinon. This chirality upsets the spin correlations, and costs energy. 
The chiral disturbance can be repaired by a second spinon nearby, which 
according to my intuition could be of either chirality. The confining force 
can hence be traced back to the appearance of chirality in the liquid as one 
pulls spinons apart. I conjecture that similar forces Dlay a central role in the 
pairing of charge carriers in CuO super conductors. 113 



11. CONCLUSION 

There are many interesting and unresolved questions regarding the 
generic spin 1 liquid state in two dimensions. The spin 1 chirality liquid 
introduced here offers a paradigm for this state and hence provides a frame- 
work to formulate, and possibly resolve, some of these. 
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